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Abstract 

Single particle electron cryomicroscopy (cryo-EM) allows for structures of proteins and protein complexes to 
be determined from images of non-crystalline specimens. Cryo-EM data analysis requires electron microscope 
images of randomly oriented ice-embedded protein particles to be rotated and translated to allow for coherent 
averaging when calculating three-dimensional (3D) structures. Rotation of 2D images is usually done with the 
assumption that the magnification of the electron microscope is the same in all directions. However, due to 
electron optical aberrations, this condition is not met with some electron microscopes when used with the settings 
necessary for cryo-EM with a direct detector device (DDD) camera. Correction of images by linear interpolation 
in real space has allowed high-resolution structures to be calculated from cryo-EM images for symmetric particles. 
Here we describe and compare a simple real space method, a simple Eourier space method, and a somewhat 
more sophisticated Eourier space method to correct images for a measured anisotropy in magnification. Eurther, 
anisotropic magnification causes contrast transfer function (CTE) parameters estimated from image power spectra 
to have an apparent systematic astigmatism. To address this problem we develop an approach to adjust CTE 
parameters measured from distorted images so that they can be used with corrected images. The efiPect of 
anisotropic magnification on CTE parameters provides a simple way of detecting magnification anisotropy in 
cryo-EM datasets. 


1 Introduction 

Anisotropic magnification in electron microscope images of 2D crystals was first described more than 30 years ago [1] 
but more recently has not been detected with most modern electron microscopes used under conditions typical for 
cryo-EM with film or charge-coupled device (CCD) cameras. Direct detector device (DDDs) cameras have revolution¬ 
ized single particle electron cryomicroscopy (cryo-EM) [2, 3]. Two DDD manufacturers, Gatan and Direct Electron, 
have produced cameras with pixel sizes between 5.0 and 6.5 pm, which is significantly smaller than is typical with 
CCD cameras. Because DDD cameras are placed below the projection chamber of the microscope, images formed 
on the DDD have an additional magnification relative to photographic film used with the same microscope. Conse¬ 
quently, electron microscopes used for high-resolution studies with a DDD may need to be set to a lower nominal 
magnification than was previously typical [4]. At these conditions, anisotropic magnification has been detected in 
several microscopes, representing a large fraction of the instruments where this issue has been investigated. The 
phenomenon was seen first in a modern electron microscope with a 300 kV FEI Titan Krios microscope used with 
a Gatan K2 Summit DDD [5] and has subsequently been detected on other FEI microscopes, including at least a 
300 kV Tecnai Polara and a 200 kV Tecnai TF20. With a 14 pm pixel size, the FEI Falcon series of DDDs does not 
require a low magnification setting and anisotropic magnification has not been detected with any microscope used 
in combination with this camera. 

The consequence of anisotropic magnification in single particle cryo-EM is that images of molecules lying in different 
orientations on the specimen grid cannot be averaged coherently, limiting the resolution that can be obtained in 3D 
maps calculated from these data. For example, with 2 % magnification anisotropy, a particle that is 300 A long 
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would appear to be 300 A long when lying in one orientation and 306 A long when lying in another orientation. This 
effect will have more severe consequences for larger particles than for smaller particles. Anisotropic magnification 
can be detected by an elliptical appearance of the powder diffraction patterns calculated from images of a variety 
of specimens including polycrystalline gold, graphite, or thallous chloride. The observed anisotropy in magnification 
exists at the low magnification settings used to acquire images, but not at the higher magnifications used for correct¬ 
ing objective lens astigmatism. The cause of this effect has been proposed to be dirt in the microscope column that 
becomes charged and acts as an additional lens; in one 120 kV microscope the anisotropy was found to change slowly 
from 1.4 % to 2.5 % over the 7 years between 1983 and 1990 (Richard Henderson, personal communication). This 
observation, as well as our own measurements described below, suggest that anisotropic magnification is stable over 
the course of weeks or months, but should likely be measured periodically with any microscope. The dependence of 
magnification anisotropy on the nominal magnification of the microscope explains why the issue was not detected in 
current microscopes when higher nominal magnifications were typical. 

In order to calculate high-resolution maps from cryo-EM images the contrast transfer function (CTF) of the mi¬ 
croscope must be corrected. Most algorithms currently in use correct for the CTF during calculation of the 3D 
map. Introduction of anisotropic magnification after correction of objective lens astigmatism causes Thon rings from 
images to be elliptical, even when there is no objective lens astigmatism present. CTF parameters determined from 
elliptical Thon rings will therefore suggest the presence of astigmatism. A consistent ellipticity for Thon rings from 
images despite several attempts at correcting astigmatism at high magnification could indicate anisotropic micro¬ 
scope magnification. Similarly, the problem may be detected by the presence of systematic astigmatism in CTF 
parameters for datasets obtained over several EM sessions where different amounts of astigmatism are expected. An 
electron optical method for removing the effect has not yet been described. Here we describe and compare algorithms 
for computationally correcting the effects of anisotropic magnification on EM images. We also develop a method 
for recovering true CTF parameters from CTF parameters calculated from images with anisotropic magnification. 
The purpose of this manuscript is to increase awareness that the potential for anisotropic magnification exists, il¬ 
lustrate how this distortion can be detected and quantified, and demonstrate how it can be corrected computationally. 


2 Methods and Results 

2.1 Measurement of magnification anisotropy 

To determine precisely the magnification of a FEI TF20 microscope operating at 200 kV and equipped with a Gatan 
K2 Summit direct detector, we acquired 99 movies containing crystalline thallous chloride particles. Anisotropic 
magnification parameters measured from these thallous chloride images can subsequently be used to correct cryo- 
EM images of biological specimens. The movies consisted of 30 frames collected at a rate of 2 frames sec“^, 5 
e“pixel”^sec“^, and 1 e” A”^ sec”^. From these movies, 558 thallous chloride particles that showed clear interatomic 
planes were selected (Fig. lA). The thallous chloride lattice has a spacing of 3.842 A. Consequently, the power 
spectrum from each crystal is expected to show a spot at distance {p • A^/3.842 A) pixels from the origin, where p 
is the pixel size in Anstroms, and N is length in pixels along each edge of the image (Fig. IB). The average of 
power spectra from images of particles is expected to produce a ring with this distance as the radius (Fig. 1C). 
However, the resulting average of power spectra had a slightly elliptical appearance. From 4096 x 4096 averaged 
power spectra of images, each containing several thallous chloride particles, we carefully recorded the lengths and 
angles of 209 vectors, d, from the origin of the pattern to the diffraction ring. Fig. ID shows a plot of the lengths 
of the vectors as a function of the angles they make with the positive /c^^-axis of the power spectrum. The plot has 
a sinusoidal appearance, indicating that the diffraction ring from thallous chloride was indeed elliptical rather than 
round. Consequently, we fit the data to the equation 

1^1 ^ |ri| + |r2| +cos(2 ■ [0^ - 6>^,]) (|ri| - [ral) 

where |d| is the distance from the origin of the diffraction pattern to a point on the elliptical diffraction pattern, ri 
and T 2 are the axes of the ellipse, Or^ is the angle between the positive /c^^-axis and ri, and Od is the angle between 
d and the /c^^-axis. The fit revealed that ri had a length that was 1.02 times the length of r 2 and is 1.3° from the 
kx-diX.is of the pattern. 
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Figure 1: Anisotropic magnification in images. A, A thallous chloride particle in a 512 pixel x 512 pixel image. Interatomic planes 
are apparent in the image. The scale bar corresponds to 100 A. B, The calculated power spectrum from the image in part A shows 
diffraction peaks. The expected distance from each peak to the centre of the image is 1.45 A • 512/3.842 A pixels. C, The average of 
558 power spectra from thallous chloride particle images. The pattern is elliptical, not round, as seen by overlay of a perfect square 
that touches the pattern at its sides but not the top or bottom. D, A plot of the distance of the pattern from the origin of the power 
spectrum versus the angle of each point from the fcx-axis in the power spectrum from the average of 99 full frame power spectra padded 
to 4096 X 4096 pixels. 


2.2 Correction of images for anistropic magnification 

Anistropic magnification can be corrected in images by appropriate stretching or contracting along ri or r 2 . Stretch¬ 
ing an image along one direction is equivalent to contracting its Fourier transform in the corresponding direction and 
vice versa. This equivalence is due to the similarity theorem [ 6 ], which states that if F{kx) is the Fourier transform 
of /(x), then is the Fourier transform of f{ax)^ where a is a constant and a 7 ^ 0. Consequently, four possible 

approaches were explored to correct images for anisotropic magnification: ( 1 ) A real space stretch along r 2 to make 
it equal to ri, (2) a real space contraction along ri to make it equal to r 2 , (3) a Fourier space contraction along r 2 
to make it equal to ri, and (4) a Fourier space stretch along ri to make it equal to r 2 . The effect of stretching or 
contracting on the position of a point in an image or Fourier transform can be expressed in matrices as a transfor¬ 
mation Fiani = ^ani^ani^Jni^ where Rjni Rani lOtate poiuts about the Origin by the angles —Oani and Oani^ 
respectively, and Sani applies a stretch or contraction along the new x-axis. These matrices are defined by: 


cos 0 dfii 

sin6>am 


sin ^ani 
COS 0dYii 


and So 77 , 7 , 


0 

1 


( 2 ) 


With these operations, the position of any point in the image p = [x^y] or Fourier transform p = [kx^ky] is 
transformed as p' = EaniP- Correcting anisotropic magnification requires applying the same transformation, but 
using a constant 1/a where the microscope caused a distortion with magnitude a. In each case, correcting anisotropic 
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magnification requires interpolation because the positions of the transformed points of the image will in general not 
fall on previously sampled image points. Real space linear interpolation is a simple interpolation scheme (Fig. 2A). 
In linear interpolation in real space in one dimension, the value of the function f{x) is estimated from the nearest 
known values of the function, /(xq) and /(xi), as f{x) ~ + xi-xo other words, the unknown 

value of the function at the point x is treated as a weighted average of the known values of the function at the points 
xq and xi. This approach can be extended two dimensions as bilinear interpolation, where f{x,y) is estimated from 
f{xo,yo), fixi,yo), fixo,yi), fixi,yi) as 


f{x,y) 


X — Xq 

Xi - Xo 
X — Xq 
Xi - Xo 


y-yo nf N , 

- j{xo,yo) + 

yi -yo 
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f{xi,yo)^ 


(3) 


Numerous more sophisticated real space interpolation schemes may also be used. Bilinear interpolation can also 
be performed in Fourier space. Interpolation in Fourier space benefits from padding images with zeros before the 
Fourier transform to increase the sampling of the Fourier transform [6]. In all of our Fourier space interpolation 
calculations we padded images to twice their original dimensions before performing the Fourier transform. Improved 
interpolation in Fourier space can be achieved with a sine interpolation approach (Fig. 2B). In Fourier space, each 
Fourier component on the regularly sampled lattice contributes to all off-lattice points as a normalized sine function 
or sin{7rAk)/{7rAk), where Ak = k — ki is the distance between the off-lattice point and the on-lattice point in the 
Fourier transform [6]. The value of the off-lattice point in a Fourier transform, F{k) is 


I 

F{k) = F{ki)sinc(TrAk) (4) 

where the summation is performed over all / on-lattice points. Summation over all Fourier components in an image 
for each interpolated point can be computationally expensive. Instead, an approximation of F{k) can be obtained 
by performing the summation over a limited set of points that are within some specified distance to the interpolated 
point. We elected to use on-lattice points that were 3 or fewer pixels away from the interpolated point as a compro¬ 
mise between accuracy and speed. 



Figure 2: Comparison of linear interpolation in real space and sine interpolation in Fourier space. A, In real space the 
value of an off-lattice point is estimated by linear interpolation, where a weighted average of the nearest on-lattice points is obtained. 
B, In Fourier space, each on-lattice point produces a normalized sine function that is 0 at all on-lattice points and non-zero at all 
off-lattice points. A point is interpolated by adding all of the contributions from all of the on-lattice points or approximated by adding 
the contributions from nearby points. The performance of this sine function interpolation is improved by padding images with zeros 
before Fourier transforming in order to increase the sampling of the Fourier transform. 
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Every interpolation approach causes different artefacts in the image, the Fourier transform of the image, or both. 
Contraction in real space (approach 2) and stretching in Fourier space (approach 4) both would cause values from 
outside of the image to be pulled into the image frame. Creating artefacts at the edge of a real space image can 
complicate the process of floating images (setting the perimeter pixels to a mean of 0) and normalization based 
on the standard deviation of the image background. Consequently, we elected not to pursue approaches 2 and 4. 
Bilinear interpolation in real space (approaches 1 and 2) causes a banding pattern in the real space images, as shown 
for an image stretched in real space (Fig. 3Ai, red arrows). This artefact arises because the interpolated value is 
always less extreme than the values from which it is interpolated (Fig. 2A). Stretching of the image in real space 
also leads to some information being shifted to higher and lower frequency in the Fourier transform of the image 
(Fig. 3Bi, green arrows). This artefact is a result of the shifting of peak positions that occurs when performing 
linear interpolation. Contraction in Fourier space (approach 4) does not lead to any notable artefacts for the result¬ 
ing image for either bilinear interpolation (Fig. 3Aii) or sine function interpolation (Fig. 3Aiii). It does, however, 
cause undefined values from outside of the Fourier transform to be pulled into the high-frequency edges of the Fourier 
transform (Fig. 3Bii and iii, blue arrows) but does not lead to any other noticeable banding artefacts in power spectra. 



Figure 3: Application of methods for correcting anisotropic magnification of images. Ai Stretching images by 2 % in real 
space causes a banding pattern artefact to appear on the image at a frequency of 2 bands per 100 pixels (red arrows). Stretching the 
image by contracting the Fourier transform with bilinear interpolation (ii) or sine interpolation (iii) causes no noticeable artefacts in 
the image. Bi, The real space stretch also causes part of the power of the diffraction peaks to be moved to higher and lower frequencies 
in the power spectrum of the image (green arrows). Contraction of the Fourier transform by either bilinear interpolation (ii) or sine 
interpolation (iii) requires that pixels be created to fill in the edge of the Fourier transform (blue arrows) but does not appear to cause 
any other notable artefacts in the power spectra from individual images. Ci, The average of power spectra from 558 uncorrected thallous 
chloride particle images shows no artefacts, other than the anisotropic magnification, ii, The average of power spectra from 558 thallous 
chloride particle images corrected for the 2 % magnification anisotropy by stretching the real space image shows several severe artefacts, 
including a non-flat background in the power spectrum (blue arrows) and movement of information to higher and lower frequency in the 
direction of the correction (green arrows), iii, The average of power spectra from 558 thallous chloride particle images corrected for the 
2 % magnification anisotropy by bilinear interpolation to contract the Fourier transform reveals a strong banding pattern (4 bands per 
100 pixels due to padding) that was not apparent from a single power spectrum (red arrows), as well as the edge pixels seen before (blue 
arrows), iv. The banding pattern is less intense when sine interpolation is used to contract the Fourier transform. 


Inspection of the average of the 558 power spectra from thallous chloride particles is a sensitive way to detect system- 
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atic artefacts in power spectra. The average of power spectra from the original images does not have any apparent 
artefacts other than the anisotropic magnification (Fig. 3Ci). The average of power spectra from images corrected 
by stretching and bilinear interpolation in real space shows severe artefacts, such as a non-uniform background 
(Fig. 3Cii, blue arrows) and shifting of power from the diffraction peaks to both higher and lower frequencies in 
the direction of the correction (Fig. 3Cii, green arrows). The average of power spectra from images corrected by 
contraction in Fourier space shows a banding pattern from the finite inclusion of Fourier terms that is more intense 
for bilinear interpolation than for sine interpolation (Fig. 3Ciii and iv, red arrows). As a result of padding images 
before Fourier transform, the frequency of these bands is twice the frequency of the artefact bands in the real space 
interpolated image. The average of power spectra from images corrected by contraction in Fourier space also shows 
the edge effects from contraction seen in a single power spectrum (Fig. 3Ciii and iv, blue arrows). Consequently, we 
concluded that the small additional computational expense of the more sophisticated Fourier space sine interpolation 
scheme is justified. Undoubtedly, it would also be possible to improve the real space interpolation approach to the 
point where it matched the performance of the Fourier space sine interpolation approach: for example by padding 
the Fourier transforms of the images to increase the sampling density or by using cubic interpolation rather than 
linear interpolation. After correction of images by contraction of Fourier transforms, the powder diffraction ring from 
thallous chloride is found at a distance of 1546 pixels from the origin in the 4096 x 4096 Fourier transform. This 
distance corresponds to a pixel size of 1.45 A and a magnification of 34,383 x on the DDD sensor. 

2.3 Effect on contrast transfer function parameters 

For high-resolution maps to be calculated from cryo-EM data, images must be corrected for the effects of the contrast 
transfer function (CTF). CTF parameters are typically measured from the power spectra of images and used to correct 
for the CTF during calculation of the 3D map. A stretch or contraction of an image causes an analogous change in 
the image power spectrum that affects the contrast transfer function (CTF) parameters determined from the power 
spectrum [7, 8]. The CTF is described by the equation: 


CTF{f) = -Wphase sin (x) - w amp cos (x) 


(5) 


where Wamp is the fractional amount of amplitude contrast (^0.07 for a cryo-EM image [9]). The parameter Wphase^ 
the fractional amount of phase contrast, is given by and 


ttA . ttA^ ^ 
X — y2 ' JF ' 


( 6 ) 


where A is the wavelength of electrons in Anstroms, / is the frequency in the Eourier transform in Anstroms”^, Az 
is the defocus of the microscope in Anstroms, and Cs is the spherical aberration of the objective lens in Anstroms. 
Eor the resolutions currently of interest in biological cryo-EM (up to ^2 A) and with the defocuses typically used 
(tens of thousands of Anstroms), the behaviour of the CTE is due almost entirely to the first term in equation 6. 
Consequently, a stretch of the image or contraction of the Eourier transform by a factor a will modify the apparent 
defocus in the image by a factor of 1/a^. This phenomenon can be observed in Eig. 4A, which compares, for a 200 kV 
microscope with a Cg of 2 mm, a ID CTE with Az =10,000 A (black line), a ID CTE with Az =10,000 A but with the 
x-axis contracted by a factor of 1.1 x (broken blue line), and a ID CTE with Az =10,000/1.1^ A= 8,264 A (red line). 


Eor a 2D image two defocus values, Azi and /S.Z 2 ^ and an angle of astigmatism, are needed to describe the CTE 
(Eig. 4B). The angle of astigmatism (j)ast is defined as the angle between the semi-axis of the Thou ring with defocus 
Azi and the /c^^-axis of the Eourier transform. The parameters Wphase^ "^amp^ A, and Cg are microscope specific while 
Azi, A 2 : 2 , and (j)ast can change with each image. Equations 5 and 6 may be used to describe a 2D CTE if one defines 
Az in equation 6 as the effective defocus in the direction <p as given by: 

^ Azi + A 2;2 + cos(2 • (0 - (l)ast)) (Azi - A 2 ; 2 ) 

During analysis of a dataset of cryo-EM images acquired on a EEI Company 200 kV EEC microscope equipped with 
a Gatan K2 Summit DDD, we plotted the measured contrast transfer function defocus parameters Azi versus Az 2 
(Eig. 4C). Note that the choice of which defocus to define as Azi and which to define as Az 2 is arbitrary and can be 
reversed by changing (past by 180°. The plot showed an apparent systematic astigmatism: for randomly introduced 
astigmatism, one would expect points in the plot to fall on the line Azi ^ Az 2 . In contrast, the points fall on a 
line that has a slope that indicated a constant astigmatism of approximately 2 %, consistent with the measured 
magnification anisotropy of the the microscope. Plotting Azi versus Az 2 and looking for deviation from the line 
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Azi = Az 2 with the form Azi ~ aAz 2 or Azi ~ presents a straightforward way of detecting anisotropic 

magnification in a microscope. However, if the anisotropic magnification factor a is small compared to the variance 
of the ratio of AzxjAz 2 due to errors in adjusting the microscope objective lens astigmatism, it will not be possible 
to detect anisotropic magnification in this way. 
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Figure 4: Effect of anisotropic magnification on CTF parameters A. The effects of stretching or contracting a Fourier transform 
on the CTF can be corrected to high precision by adjusting the defocus parameter. B, Contrast transfer function (CTF) parameters 
that are unique for each image are the two defocus parameters A 2 ;i, A 2 ; 2 , and the angle of astigmatism <^ast- C, A plot of Azi vs I\z 2 
for CTF parameters determined from images that suffer from magnification anisotropy show an apparent systematic astigmatism, with 
Azi consistently being larger (or smaller) than A^2- Note that the ratio of Azi to can be reversed by redefining which defocus 
corresponds to which parameter and adding ^ to bast- C, The apparent systematic astigmatism induced by anisotropic magnification 
can be removed by correcting CTF parameters as described in the text. 


Correction of images for anisotropic magnification before measuring CTF parameters removes this systematic devi¬ 
ation from the Az\ = Az 2 line. However, it is also possible to correct CTF parameters that have been measured 
from images with anisotropic magnification. It is advantageous to measure CTF parameters from images with mag¬ 
nification anisotropy and then correct the parameters for three reasons. First, computationally correcting whole 
images and movies, rather than individual particle images, could be computationally expensive. Second, storing 
computationally corrected whole images and movies places significant demand on data storage resources, which may 
already be strained by DDD data collection. Finally, correcting images for anisotropic magnification can introduce 
image artefacts that may decrease the accuracy with which CTF parameters can be calculated. Recovery of true 
parameters from images with a known amount of magnification anisotropy relies on representing the CTF as a dis¬ 
tortion of the unit circle in Fourier space, -h = 1, to an ellipse (Fig. 4B). The unit circle may be represented in 
matrix notation as 


k'^k = 1 


(8) 
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where = [kx ky and k 


kii 


. As seen in Fig. 4B, an ellipse can be used to represent the CTF parameters 


Azi, Az 2 , and (past- Defining k = 'Ectf^ where Ec 


CTF = ^CTF^CTF^CTF 




with 


R. 


CTF — 


COS 00,5^ 


and S 


CTF — 


0 




and substituting into equation 8 gives the equation of the ellipse as 


^'EJ.tf^ctf^ = 1. (9) 

^CTF describes a rotation by the angle of astigmatism past from the /c^^-axis, Sctf describes a stretch or compression 
along the kx- and ky-diX.es by Azi and Az 2 , and Hctf rotates the ellipse axes back to their correct angles. The 
magnitudes of the two defocus values {Azi and AZ 2 ) can be recovered by determining the two eigenvectors of 
'E^j.p'EcTF- The angle of astigmatism past can be calculated as the angle that the eigenvector representing Azi 
makes with the kx-Sixis. By definition, the eigenvectors of a matrix Y obey the equation Yx = Ax. That is, upon 
multiplication by Y the eigenvectors change only their magnitude, not their direction. Eigenvectors of a 2 x 2 matrix 
can be conveniently obtained using algorithms in standard numerical analysis libraries such as LAPACK, NumPy, 
or in MATLAB. As described earlier, magnification anisotropy further distorts the ellipse by the transformation 
E^ni = RaniSaniRani where Rani describes a rotation by the angle of magnification anisotropy Oani from the kx- 
axis, and Sani describes a stretch or compression along the kx-dxis by the magnitude of magnification anisotropy a, 
where Rani and Sani are the shown in 2. Defining i = Eaniiai and substituting this definition into equation 9 gives 
the equation of the new ellipse as 


m^EVE^r^EcTFEanim = 1. (10) 

The eigenvectors of Ej^^E^^^ECT^^Eani correspond to the apparent CTF parameter values, Az[^ and p^st^ 

measured from the power spectra of images where there is aniostropic magnification. From these values and values 
of a and Oani measured from a powder diffraction pattern (see above) it is possible to calculate the true values of 
Azi, Az 2 , and past from the eigenvectors of 

E T TTt _ 17 >—TfT'/T T 7 ^/ IT' —1 /"I 1 ^ 

CTF^CTF — ^ani^CTF^CTF^ani 1 -'--'-) 

where ctf''R-ctf' and EV = The matrices Rctf', and Sctf' are the 

same as RctFj R-ctF’ and Sctf, but use the apparent CTF parameters measured from distorted images rather 
than the true CTF parameters, and 


0-1 ^ [l/n 0 

^ani Q ^ 

This approach for correcting CTF parameters measured from images with magnification anisotropy was implemented 
in a standalone program that operates on Relion .star files [10]. As can be seen in Fig. 4D, the method brings points 
in a plot of Azi versus Az 2 back to the Azi = Az 2 line. The method is equivalent to correcting an image for 
magnification anisotropy and subsequently measuring CTF parameters from its power spectrum. 


3 Discussion 

As seen in Fig. 4C, a plot of Az\ versus Az 2 from a dataset of micrographs can reveal the presence of anisotropic 
magnification in the microscope. If all of the images in a dataset were obtained from a single EM session where 
the objective lens stigmator was adjusted once, it would not be unusual to find the points of this plot fall off the 
line Az\ = Az 2 . However, if the deviation of points from the line is due to anisotropic magnification, the positions 
of the points should have the form Az\ ^ aAz 2 or Azi ^ yAz 2 . Deviation due a constant astigmatism should 
result in point positions with the form Azi ^ Az 2 + where x is the the amount of astigmatism. Once anisotropic 
magnification has been detected, the amount and extent of anisotropy can be measured precisely using a diffraction 
standard such as the thallous chloride shown above. The interpolation scheme proposed here is not entirely free 
of artefacts. Consequently, it would be better to correct anisotropic magnification with improved electron optics, 
rather than correct the effects of anisotropic magnification computationally. If correction by interpolation must be 
performed, the best approach would be to incorporate the correction into the 3D map calculation software in order 
to avoid performing interpolation twice with the same images. 
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4 Software availability 

All of the original software described above is available from https://sites.google.com/site/rubinsteingroup/direct- 
detector-distortion 
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